EFFECTIVE FIELD THEORY FOR NUCLEAR PHYSICS 



DAVID B. KAPLAN 

Institute for Nuclear Theory, University of Washington, Seattle, WA, 98195-1550, 

U. S. A. 

I summarize the motivation for the effective field theory approach to nuclear 
physics, and some of its recent accomplishments. 



1 Why effective field theory, why nuclear physics? 

Low energy data is generally insensitive to the details of interactions at short 
distance. It is therefore difficult to learn about short range interactions; yet 
by the same token, complete knowledge of the physics at short distances is 
not required for an accurate understanding of experiments. Effective field 
theory exploits this fact. The effects of nonlocal interactions at short distance 
may be represented in terms of local operators in a derivative expansion — 
the effective Lagrangian. The higher an operator's dimension, the smaller 
the effect it has on low energy physics, and hence one can obtain a useful 
phcnomenological theory by retaining operators only up to some dimension, 
fitting their coefficients to data. Some effective theories are quite useful, such 
as chiral perturbation theory; some are wildly successful, such as the standard 
model of particle physics. In this talk I will discuss a new application currently 
being developed, nuclear effective theory. 

The utility of effective field theory (EFT) depends on the existence of an 
energy gap so that "short" and "long" distance physics can be distinguished. 
It is probably not a useful technique for describing turbulence, or protein fold- 
ing, for example. In 1990, Weinberg suggested that nuclear physics could be 
a subject that would benefit from an EFT treatment In nucleon-nucleon 
interactions, one can identify the low scales to be m T = 140 MeV, and the 
nucleon momentum (pp ~ 280 MeV in nuclear matter), while the high scales 
would be the masses of the vector mesons and higher resonances. The A con- 
tributes to NN scattering at a momentum scale ~ y/2M^{M/^ — Mn) = 750 
MeV, comparable to the vector mesons. While Weinberg's original proposal 
has been shown to not be consistent, a new approach-has been developed over 
the past nine months that appears very promising □ (see also Ref. a). It is 
this approach that I will focus in this talk. 

What are the benefits of EFT for nuclear physics? The standard starting 
point of nuclear theory for the past half century has been to take the NN 
scattering phase shifts as input, and then construct a Schrodinger potential 
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that can reproduce those phase shifts. Using this potential on can solve the 
Schrodinger equation for few body systems, or use various many-body ap- 
proximation schemes (Bruckner theory, the shell model, etc.) to solve larger 
nuclei. The main weakness of this approach is that there is no systematic way 
to improve results, or to have confidence in their level of accuracy. The lack 
of a systematic expansion becomes evident in many-body calculations when 
one finds dependence on "off-shell amplitudes" . As is well known, off-shell 
amplitudes are arbitrary and can be changed at will by field redefinitions. 
That a purportedly physical quantity depends on the off-shell definition of 
an amplitude is a sure sign that a consistent expansion is lacking. As I will 
discuss below, a systematic expansion is obtainable in the EFT approach to 
nuclear physics, so that one never need to specify off-shell matrix elements. 

EFT offers other, practical, advantages as well: by replacing complicated 
models of short distance physics by separable contact interactions (in a justi- 
fiable fashion), calculations are vastly simplified. In the two nucleon sector, 
many computations have been performed analytically, while in the three nu- 
cleon system, the Fadeev equations are reduced to a one dimensional integral 
equation quickly solvable on a desktop computer. Furthermore, implemen- 
tation of relativity as well as chiral and gauge symmetries is simpler using a 
local effective Lagrangian, compared to a more traditional nucleon potential 
approach. 

2 A systematic expansion 

Developing a systematic expansion for an EFT for nuclear physics is less 
obvious than for more familiar applications. Typically one introduces all op- 
erators allowed by low energy symmetries, introduces couplings which depend 
inversely on the high energy scale A to the power appropriate for the dimen- 
sion of the operator, and then compute the Feynman amplitude to the desired 
power of p/A or m/A, where p is a external momentum and m is a light par- 
ticle mass. If performed sensibly, the computation never introduces powers of 
A in the numerator and so at which order an operator enters depends solely 
upon its dimension. 

The lowest dimension contribution to NN scattering at low energies would 
come from the operator Cq(N^ N) 2 (ignoring spin and isospin indices) where N 
is the nonrelativistic nucleon doublet field and Co is a coupling constant with 
dimension (mass) -2 . If EFT for nuclear theory worked in the accustomed way, 
then one would expect Co ~ l/(MA)_where A is the vector meson mass scale, 
and M is the nucleon mass (see Ref . □ for a discussion of the M dependence) . 
One would then insert this Co perturbatively: each insertion would contribute 
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an additional power of MpCo ~ p/ A to the amplitude, which is small at low 
energy. 

The reason why this cannot be the whole story behind NN scattering is 
evident if one looks at the NN phase shifts. In the 1 5'o channel np scattering 
exhibits a scattering length a = — 1/(8 MeV), while in the 3 Si channel the 
np scattering length is a = 1 / (40 MeV) . These length scales are obviously 
neither the inverse pion or vector meson masses, but must arise through non- 
perturbative dynamics. In the potential model approach, they are due to a 
near cancellation between potential and kinetic energy. The effective field 
theory description instead resembles a condensed matter system near a phase 
transition, a system that similarly exhibits length scales (correlation lengths) 
much longer than the fundamental scale of the system. In other words, one 
will have a quantum field theory tuned to lie near a nontrivial fixed point, 
where the fixed point corresponds to infinite scattering length. The EFT is 
then necessarily nonperturbative, and the actual dimension of operators can 
be different than their naive dimension, with power law corrections to the 
scaling of coupling constants. 

To exhibit this fixed point behavior, consider NN scattering at momenta 
far below m„. The scattering amplitude is then well described solely in terms 
of the scattering length a, 

47T 1 , . , 

M 1/a + ip 

This can be reproduced by the effective theory consisting of nonrelativistic 
nucleons interacting via the Co(iV'iV) 2 term. By iterating bubble diagrams 
with a Co at each vertex, one finds 

where I(p) is the linearly divergent loop integral 

f d 3 q M 

J [2-k) 6 p z - q z 

If one evaluates this diagram using dimensional regularization and the PDS 
subtraction scheme introduced in Rcf. larj-similarly, if one performs a mo- 
mentum space subtraction at p 2 — ~/i 2 Jiltffl) one finds that in terms of the 
renormalized coupling Cq(ja), the amplitude is 



•Aeff — Ti/T 4-tt , , • W 

M TRMrf + M + l P 
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Comparison with eq. implies 

Note that the dimensionless coupling constant defined as Cq = — M ^,Cq / (4w) , 
satisfies the renormalization group equation 

^ = 0,(1 - C ) . (6) 

This equation features a nontrivial UV fixed point at C'o = 1, corresponding to 
1/a = 0. At this fixed point Coin) oc 1///, so that in fact Co(N'N) 2 behaves 
as a marginal operator, instead of irrelevant. 

We can now do some reverse engineering and identify a consistent power 
counting scheme: Since Co(/i) oc 1//Lt, and each loop brought with it a factor 
of /i or p, we can define an expansion variable Q with fi ~ p ~ Q. Then 
we see that every term in the bubble sum we just performed was order 
This is our leading order amplitude. The final expression Eq. (|l|) implies that 
1/a ~ Q\ our expansion is only valid for systems with large scattering length. 

To proceed to higher order we consider higher derivative contact interac- 
tions 

Co + C 2 p 2 + C 4P 4 + ... (7) 

as well as interactions involving pions and photons. Simple power counting 
rules emerge, allowing one to systematically expand the amplitude (the KSW 
expansion): 

1- C 2nP 2n ~ Q"- 1 ; 

2. - Q; 

3. V ~ Q , dt ~ Q 2 for derivative couplings of the pion or photon; 

4. Loop integrations d 4 q count as ~ Q 5 if go ~ Q 2 /M, otherwise as Q 4 ; 

5. propagators count as Q~ 2 . 

Armed with these rules, one computes physical observables to a given order 
in Q. By working consistently to a given order, one is ensured of always 



"Rule (1.) applies to s-wave scattering; forr* more detailed discussion of power counting 
including arbitrary partial waves, see Ref. 0. The second option in rule (4.) occurs at 
NNLO where radiation pions first appear. 
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being able to renormalize the theory (one has enough couplings to absorb 
divergences) and no observable will depend on off-shell matrix elements. 

The above rules show that effective nuclear theory contains chiral pertur- 
bation theory with its usual perturbative treatment of the pion...one never 
needs to sum up pion exchange to all orders. That is useful in practice, as 
it allows for analytical calculations. It is also satisfying theoretically, since 
it would be hard to imagine chiral perturbation theory to be relevant if one 
needed to treat pions nonperturbatively. 



3 Application to a toy model 

Before proceeding to the real world, it is worthwhile to analyze a toy model 
that exhibits two distinct length scales and which can be solved analytically. 
That toy model is a pair of nucleons interacting via delta shell potentials 

Tn Til 

V(r) = g p ^8{r - l/m p ) + g^S(r - l/m«) . (8) 

The couplings g p , g n are dimensionless and normalized so that g = — 1 is a 
strong coupling just capable of producing a bound state. In order to have 
a power counting similar to what is found in the real world, I define g^ = 
— m^/A, where A has dimensions of mass. 

The phase shift 6 may be computed analytically in this model, and can 
be taken to be a function of the five variables {m p , A, to w , a,p}, where a is the 
scattering length and p = \JME cm . If we perform the rescaling 

— > em^ , p — > ep a — > —a , (9) 

and expand the phase shift in powers of e as 

S(m p ,A, errin, a/e, ep) = S Q + 8\e + ^ 2 ( 2 + ■ ■ ■ (10) 

then, following the discussion in the previous section, we would expect the 
functions 5 n (m p , A, m„ a,p) to correspond to the n th term in the KSW ex- 
pansion of the effective theory. In fact, this is correct. When the short range 
part of the interaction in Eq. (^J) is replaced by contact interactions, 

V e ff = C + C 2 p 2 + C^p 4 + . . . + 9tt-^K t ~ V"V) , (11) 

and one expands the amplitude using the rules of the previous section, one 
finds that appropriate choices of the C2„ couplings allow one to recover the 
expansion Eq. ([h]) of the exact result, order by order []. 

6 The procedure involves an expansion of Co as discussed in Refs. fJB 
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Figure 1. Leading (dotted)and subleading (dashed) EFT result for the toy model, compared 
with the exact answer (solid), plotted as degrees vs. p = \/ME cm in GeV. 
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Figure 2. Fourth order EFT result (dashed) compared with the exact result (solid), plotted 
as degrees vs. p = \/ME cm in GeV. 



How quickly does the series converge to the true phase shift? I take 
m„ = 140 MeV, m p — 770 MeV. In order to have the toy model resemble 
realistic np scattering in the channel, I choose a = — 23 fm, and set 
A = 500 MeV so that the "pion" interaction in the toy model resembles the 
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Figure 3. The effective range approximation for the phase shift in the toy model (dotted) 
compared with the exact result (solid), plotted as degrees vs. p = \/ME cm in GeV. The 
approximation breaks down at p ~ . 



true one-pion exchange potential [j. In Fig. |l] I have plotted the first two 
terms in the KSW expansion (dashed lines) against the exact result (solid 
line). The lowest order (LO) result gets the scattering length exactly, but 
quickly fails to describe the "data". The next-to- leading order (NLO) result 
is a vast improvement, working crudely up to p ~ m p . In Fig. |2| I show the 
result at fourth order in the KSW expansion; it exhibits excellent agreement 
up to p ~ m p , and shows not only that the expansion is converging rapidly 
(in spite of the enormous difference between LO and NLO results), but also 
that the radius of convergence in momentum is set by m p and not m w . For 
comparison, I have plotted the effective range result in Fig. ||. It does an 
excellent job at describing the true phase shift for p < m x (better than the 
NLO EFT result), but fails as one would expect above p ~ 

Recently several papers have appeared suggesting that pions must be 
included, nonperturbatively in order to obtain convergent answers above p ~ 
nix B'EhIJ. This does not appear to be the case in the toy model, and I don't 
believe it is true for the real world. Instead, I think one gets into trouble 

c By choosing A = 500 MeV, we get g n = — mvr/A = —0.28, which is 28% of the critical 
value that would give rise to a boundstate due to pion exchange alone; the real one-pion 
exchange Yukawa potential in the 1 Sq channel has a strength which is similarly 28% of the 
critical value. Furthermore, for both the toy and realistic pion potentials, the maximum 
phase shift in the absence of additional interactions is ~ 12° . 
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Figure 4. NN phase shifts from Ref. □. The dot-dash and dotted lines are the LO and NLO 
EFT results respectively. The solid lines are from the Nijmegen partial wave (multi-energy) 
analysis; the data points for ei are from a single-energy analysis H3- 



by trying to exactly reproduce the effective range expansion parameters a 
and r at finite order in the KSW expansion. That is because one is forcing 
the short range physics to account for the contributions to these observables 
from neglected higher order pion exchange. By incorrectly determining the 
short distance physics (which is finely tuned), one runs into trouble at higher 
momentum. In the toy model considered here, choosing the CVn coupling to 
exactly reproduce the (n + 1) term in the effective range expansion yields an 
EFT prediction that fails (at any order in the expansion) at p ~ m,. These 
issues are to be discussed at greater length in Ref. LLJ. 

4 The real world 

I will now briefly describe the status of calculations for real two and three im- 
cleon physics using the KSW expansion. Phase shifts have been fit to NLO □ in 
both the 1 S'o and 3 Si — 3 Z?i channels, as shown in Fig. ||. Using the C param- 
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eters derived from the NN phase shifts, a number of two-nucleon processes 
have been computed. In Figs. ((5||) I show results for the electromagnetic 
form factors of the deuteron from Ref. 13. These figures give evidence that 
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Figure 5. Electric form factor of the deuteron. The dotted and solid curves are the LO 
and NLO effective field theory results respectively; th« dashed curve is the effective range 
theory prediction. The data are from various sources Q. 

the KSW expansion is converging. However, it also shows that conventional 
effective range theory gives a somewhat better fit to the data than does the 
EFT calculation at NLO. Presumably this is because the effective range ap- 
proximation is fit to reproduce low energy scattering data exactly, while the 
EFT calculation yields low energy amplitudes in an expansion in powers of 
~m w /A~l/3. 

In Fig. [?] is displayed the result from Ref. Ii3 for Compton scattering 
off the deuteron. Recently np capture has been examined, both the._Darity 
conserving part to NLO E_a, and the parity violating part to LO Ref. Ej. All 
of these calculations have been performed analytically. To give you an idea 
what is involved, I display in Fig. [| the graphs that had to be computed 
for the calculation of parity violation in np capture. Other two-nucleon EEI 
calculations usingihe KSW expansion include the deuteron anapole momcntliS 



The most interesting challenge in nuclear effective theory presently is to 
extend the analysis to nuclear matter or finite nuclei. There are a number 
of issues that arise in the three-body system that must be understood be- 
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Figure 6. Magnetic form factor of the deuteron. The dotted and solid curves are the LO 
and NLO effectivo-ifield theory results respectively; the dashed curve is the effective range 
theory prediction Q. 

fore progress can be made. The simplest three-nucleon process to consider 
is A^cLscattcring in the J — 3/2 channel. This was recently computed in 



made in understanding the more interesting J = 1/2 channel from the EFT 
perspectiveEx The problem is that the graphs one must compute must be rep- 
resented as an integral equation, and no systematic power counting has yet 
been developed , although many interesting features have been understood. 
Understanding this system promises to be rewarding. 

5 Discussion 

How does the EFT approach measure up to conventional nuclear physics, 
where one solves the Schrodinger equation with a phenomcnologically derived 
potential? So far, the LO and NLO calculations performed in the past nine 
months are not of comparable accuracy. Potential models obviously do much 
better than EFT at fitting NN phase shifts, as they typically have about 
two dozen parameters expressly tuned to fit the data. What about other 
quantities? 

• Deuteron electric moments: State-of-the-art potential models also do 
quite well at explaining quantities particularly sensitive to long distance 
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Great progress has also been 





Figure 7. Angular distribution for Compton scattering off the deuteron, for two different 
energies, fn each case the dashed (red) curves represent the LO result; the solid (blue) 
curves are the full NLO result. The dotted (fed) curves are the NLO result when nucleon 
polarizability graphs are omitted. From Ref. E3. 



physics. This includes the deuteron charge radius and magnetic moments, 
accurate to ~ 1% E3. In contrast, the EFT predictions foiLthe deuteron 
charge radius are off by 28% at LO, and by 11% at NLO 0. The reason 
for this rather poor agreement is that even low momentum properties 
get corrections order by order in the KSW expansion. The deuteron 
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Figure 8. rfleynman diagrams contributing to parity violating np capture to lowest order, 
from Ref. E3. The square vertex arises from the weak interactions. 



quadrupole moment is off by 6% in potential model calculation^; only 
the LO result from EFT has been obtained to date, and it is off by 39%. 

• Deuteron polarizability: The polarizability of the deuteron has not been 
measured, but different potential models agree with.. each other on the 
value of the scalar electric polarizability to ~ 0.3% EJ; the EFT calcula- 
tion differs from the potential model results by ~ 39% at LO, and 6% at 
NLO. 
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Figure 9. N — d scattering phase shift, J = 3/2 channel, computed in EFjPriThe only input 
are the two-body J = 1 scattering length and effective range. From Ref. 



• Deuteron Compton scattering: Compton scattering in the energy range 
shown in Fig. [?] is less sensitive to the very long range features of the 
deuteron wavelength. Predictions from potential models E3 agree with 
the data at the ~ 10% level, which is comparable to the agreement found 
in the NLO EFT calculation. 

My conclusion is that the operative expansion parameter in the KSW expan- 
sion is ~ 1/3, with LO results agreeing with data to the 30% level and NLO 
agreeing to the 10% level. NLO calculations are typically not very arduous, 
and can usually be obtained in closed, analytic form. An important NLO 
calculation yet to be performed is the quadrupole moment of the deuteron. I 
expect that NNLO calculations, expected to be accurate to ~ 3%, will prove 
to be competitive with — and in some cases, superior to — potential model 
calculations for a wide variety of two nucleon processes. Such calculations 
involve relativistic corrections and propagating pions, and have not yet been 
fully understood, although several groups are working on the problem. NNLO 
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results most competitive with potential model results will be those relatively 
insensitive to the long wavelength tail of the deuteron wave function. Perhaps 
a better understanding of how in practice to match the KSW expansion to the 
effective range expansion at very low energies could improve the predictions 
for the latter observables. 

Extending the EFT approach to few- and many-body systems remains an 
intriguing challenge, one that must be surmounted before EFT can really be of 
much use in nuclear theory. Like civilizations, physics theories progress from 
barbarism, to civilization to decadence (occasionally skipping the civilized 
state). Nuclear effective theory is still in the barbaric stage. Nevertheless, 
great progress has been made over the past nine months, and I am optimistic 
that the techniques will catch on and eventually prove themselves invaluable^. 
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